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INTRODUCTION 


Each  year  the  U.S.  Navy  is  faced  with  the  problem  of  allocating 
its  procurement  budget  to  various  items.  For  the  Navy,  the  two  primary 
items  are  ships  and  aircraft.  This  paper  will  restrict  its  attention 
to  these  two  items  and  the  procurement  budget  associated  with  them. 
Thus,  each  year  a  certain  proportion  of  that  budget  must  be  allotted 
to  the  purchase  of  aircraft,  with  the  remainder  being  spent  on  ships. 
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But  what  should  that  proportion  be?  The  answer  to  that  question 
depends  on  what  the  Navy  considers  to  be  an  optimal  force  mix  between 
aircraft  and  ships  (i.e.  ,  the  ratio  of  the  number  of  aircraft  to  the 
number  of  ships).  The  objective  of  any  budget  allocation  plan  should 
be  to  achieve  this  optimal  force  mix  ratio  between  aircraft  and  ships. 

Constraining  the  Navy’s  attainment  of  this  optimal  force  mix, 
however,  is  the  procurement  budget  available  and  the  respective  costs 
of  aircraft  and  ships  over  a  number  of  years.  If  the  available  budget 
is  small,  the  Navy  may  have  to  make  a  trade-off  between  aircraft  and 
ship  procurement,  allowing  one  force  (either  aircraft  or  ships)  to 
shrink  in  order  to  expand  the  other  and  so  obtain  the  desired  force 
mix  ratio  in  that  way.  On  the  other  hand,  if  available  budgets  are 
large,  thus  allowing  an  increase  in  both  forces,  the  Navy  could  obtain 
its  desired  force  mix  by  expanding  both  forces  by  some  ratio  and  attain 
the  goal  with  both  forces  at  a  higher  level.  In  either  case  there 
arises  the  need  to  allocate  the  available  budget  in  such  a  way  as  to 
stay  as  close  to  the  desired  ratio  as  possible  over  some  planning 
period.  While  the  Navy  does  not  specifically  plan  forces  based  on  a 
set  ship  to  aircraft  ratio,  the  historical  facts  indicate  that  a 
relatively  constant  ship  to  aircraft  ratio  results. 

Thus,  the  objective  of  the  optimal  budgeting  plan  considered 
here  will  be,  given  information  on  budgets  and  costs  over  a  planning 
period,  the  minimization  of  the  deviation  of  the  actual  force  mix 
ratio  from  the  desired  force  mix  ratio  for  all  years  in  the  planning 
period. 
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PARAMETERS 

In  developing  an  optimal  budget  allocation  many ' parameters  will  need  to  be 
taken  into  account.  An  explanation  of  several  of  the  most  important  ones  follows. 

First,  there  are  the  respective  retirement  rates  for  aircraft  and  ships. 

Ships  in  the  fleet  last,  on  average,  30  years.  So,  assuming  the  ages  of  the 
fleets  ships  to  be  spread  somewhat  uniformly  between  0  and  30  years,  one  would 
expect  to  xose  l/30  (3*3/^)  of  the  fleet  each  year.  Currently  there  are  approx— 
imately  35 0  ships  in  the  fleet.  In  order  keep  the  fleet  at  Its  current  level, 

18  (550/30)  new  ships  must  be  added  each  year.  As  for  aircraft,  the  average 
aircraft  lasts  15  years  and  so  roughly  l/l5  or  6.? %  of  the  Navy’s  air  fleet  will 
retire  each  year.  Currently  the  are  3500  aircraft  supported  by  the  Navy,  so, 
to  maintain  current  air  force  levels,  approximately  235  new  aircraft  will  be 
needed  each  year. 

Second,  there  is  the  problem  of  spiraling  aircraft  costs.  Since  aircraft 
costs  are  increasing  at  a  rate  greater  than  either  ship  costs  or  the  procurement 
budget  growth  rate ,  it  becomes  increasingly  difficult ,  in  later  years ,  to  main¬ 
tain  the  desired  force  mix  without  allowing  the  number  of  ships  in  the  fleet 
to  decline . 

A  third  parameter .bearing  on  an  optimal  budget  allocation,  is  the  leadtime 
for  delivery  of  aircraft  and  ships.  It  is  estimated  that  the  leadtime  for  air¬ 
craft  is  2  years  and  for  ships  5  years.  Hence,  a  budget  allocation  in  a  given 
year  affects  the  ratio  of  aircraft  to  ships  both  two  and  five  years  from  the 
allocation  year .  These  differential  leadtimes  lead  to  an  extensive  inter¬ 
relationship  between  the  various  years  of  budget  allocation.  For  instance, 
the  force  mix  ratio  in  year  i  is  affected  by  the  ratio  in  year  i-1  and  the 
budget  allocations  in  years  i-2  and  i-5  which  are  in  turn  affected  by  actions 
or  conditions  in  other  years  of  the  planning  period.  Thus,  looking  at  any  one 
year's  budget  and  analyzing  a  certain  allocation's  effects  on  any  one  other 
year  would  not  be  valid,  The  budget  allocation  scheme  must  be  viewed 
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system  spread  over  the  many  years  of  the  planning  period. 

Finally,  there  is,  to  reiterate,  the  procurement  budget's  affect  on  the 
desired  force  mix.  If  the  budgets  are  large,  relative  to  costs,  force  levels 
will  grow  while  if  budgets  are  small  levels  will  fall.  This  fluctuation  of 
force  levels  with  respect  to  the  available  budget  will  necessarily  need  to  be 
taken  into  account  when  determining  the  optimal  budget  allocation  (i.e.,tne 
budget  allocations  should  be  different  for  growing  and  contracting  foree  levels). 

THE  MODEL 

Now  that  some  of  the  parameters  affecting  the  optimal  budget  allocation 
have  been  explained,  a  description  of  the  model  to  calculate  the  optimal  pro¬ 
portion  of  the  budget  to  be  spent  on  aircraft  each  year,  with  remainder  being 
spent  on  ships,  can  be  made.  The  assumptions  of  the  developed  model  are: 

1)  Retirement  rates  for  aircraft  and  ships  are  constant  (i.e.,  a  fixed 
percentage  of  the  total  fleet  of  aircraft  and  ships  will  be  retired 
each  year) ; 

2)  Leadtimes  for  aircraft  and  ships  are  2  and  5  years  respectively; 

3)  There  is  a  fixed  planning  period  (some  #  of  years)  over  which  we  are 
concerned  about  the  force  mix  ratio; 

4)  The  available  budgets  and  costs  of  aircraft  and  ships  are  known  for 
all  years  in  the  planning  period. 

The  objective  in  choosing  the  proportion  of  the  budget  to  be  spent  on  aircraft , 
with  the  remainder  being  spent  on  ships,  will  be  to  minimize  the  sum  of  the 
squared  deviations  of  the  actual  force  mix  ratio  from  the  desired  force  mix 
ratio.  In  other  words,  we  want  a  series  of  proportions,  r,  ,  one  for  each- 
year  in  the  planning  period,  such  that  these  r^  minimize  the  sum  of  the  squared 
deviations  of  the  actual  ratio  from  the  desired  one.  The  reason  for  the  squared 
deviation  is  the  fact  that  ratios  either  above  (yielding  a  positive  deviation) 
or  below  (yielding  a  negative  deviation)  are  considered  equally  costly,  so  a 
deviation  of  -1  or  +1  should  be  treated  the  same .  This  is  done  through  squaring 
the  deviation.  This  method  will  force  the  actual  ratio  to  stay  as  close  as 
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possible  to  the  desired  force  mix  ratio.  This  objective  is  then  subject  to 
the  available  budgets,  costs  of  aircraft  and  ships,  .and  further  to  the  con¬ 
straint  that  the  r^  must,  logically,  be  between  0  and  1. 

Formulated  in  this  manner,  the  problem  of  finding  the  proportion  of  the 
budget  to  be  spent  on  aircraft  each  year  is  a  nonlinear  programming  problem 
with  the  r^  proportions  as  the  decision  variables.  From  here,  using  the  elements 
of  multivariable  calculus  and  optimization  theory  it  is  possible  to  develop 
a  concise  format  for  the  stated  problem, implement  it  on  a  computer,  and  solve  it 
for  any  set  of  input  data. 

INPUTS 

The  inputs  to  the  developed  model  are  as  follows: 

1)  Number  of  years  in  the  planning  period  (up  to  30) 

2)  Desired  force  mix  ratio  (aircraft  to  ships) 

3)  Current  numbers  of  aircraft  and  ships 

4)  Available  procurement  budget  (for  aircraft  and  ships)  for  each  year 
in  the  planning  period 

3)  Aircraft  cost  in  each  year  (average) 

6)  Ship  cost  in  each  year  (average) 

7)  Retirement  rates  for  aircraft  and  ships 

8)  Number  of  aircraft  ordered  1  and  2  years  ago 

9)  Number  of  ships  ordered  1,2,3(4  and  3  years  ago 

It  should  be  noted  that  inputs  8  and  9  represent  backorders  of  aircraft  and 
ships  which  will  arrive  diming  the  planning  period  but  which  are  hot  under  the 
control  of  the  optimal  budget  allocation  scheme . 

AN  EXAMPLE 

With  these  concepts  in  mind  and  with  the  developed  software,  it  is  now 
possible  to  construct  a  numerical  example  and  show  how  the  optimal  budget 
allocation  scheme  might  better  achieve  the  desired  force  mix  ratio  than  a  less 
flexible  approach. 
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Consider  this: 

1)  A  planning  period  of  20  years 

2)  A  desired  ratio  of  7  aircraft  to  1  ship 

3)  Current  force  levels  of  3500  aircraft  and  550  ships 

4)  Available  budget  starting  at  23  billion,  increasing  at  3%  per  year 

5)  Aircraft  cost  starting  at  40  million,  increasing  at  4^  per  year 

6)  Ship  cost  starting  at  552  million,  increasing  at  3%  per  year 

7)  Retirement  rates  for  aircraft  and  ships,  &.7%  and  3  >3%  respectively 

8)  Backorders  of  aircraft  and  ships  to  just  cover  retirements,  235  aircraft 
per  year,  18  ships 

The  results  of  the  optimal  budget  allocation  program  are  given  below  and  on  the 
following  page.  The  r(i)  column  represents  the  proportion  of  the  budget  which 
should  be  spent  on  aircraft  each  year.  Notice  how  close  the  optimal  allocation 
scheme  allows  the  actual  ratio  to  stay  to  the  desired  ratio.  Also  notice  the 
variablility  in  the  proportion  spent  on  aircraft.  This  is  primarily  due  to  the 
changing  ratio  of  ship  cost  to  aircraft  cost. 
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In  contrast  to  that,  consider  the  results  of  setting  r(i)  equal  to 
the  average  proportion  in  the  optimal  plan  for  all  years  in  the  planning 
period.  This  results,  see  below,  in  a  ratio  of  aircraft  to  ships  which 
varies  substantially  from  the  desired  ratio .  Thus ,  in  this  example ,  the 
efficiencies  gained  through  the  mathematical  programming  formulation  of 
the  problem  can  be  seen  through  its  ability  to  remain  very  close  to  its 
objective  while  other  methods  cannot  seem  to  equal  this  ability. 


AVERAGE  OPTIMAL  ALLOCATION  RESULTS 


Desired  Ratio  “  7.00 


Fee  tire me n  t  I < 5 to  for  Aircraft  -  0 . 06  7 
ftetj  remer.t  Rate  for  Ship's  :=  0.033 
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Though  the  average  proportion  results  presented  in  the  foregoing  example  were 
not  really  "bad,"  the  fact  is  that  even  these  results  could  not  have  been  ob¬ 
tained  without  the  optimal  budget  allocation  software.  Hence,  at  least  to  a 
certain  extent, the  preceding  example  illustrates  the  efficacy  of  the  method. 

ANOTHER  EXAMPLE 

With  the  developed  model  an  analyst  or  decision  maker  could,  with  a  good 
deal  of  accuracy,  asses  the  force  mix  impact  of  various  changes  to  the  model. 

For  instance,  the  effects  of  very  small  budgets  on  absolute  force  levels  could 
be  studied  given  that  the  force  mix  was  to  be  kept  constant.  Alternatively, 
one  could  examine  the  effects  on  absolute  force  levels  of  changes  in  the  desired 
force  mix  ratio.  As  an  example  consider  this: 

1)  The  Navy  wants  to  keep  at  least  a  500  ship  fleet 

2)  Budgets  available  are  small  (starting  a  $15  billion,  increasing  at  J %  per  yr) 

3)  The  Navy  would  like  to  maintain  a  7  to  1  ratio  between  aircraft  and  ships 

4)  Other  model  inputs  are  the  same  as  in  the  last  example 

Are  the  above  criteria  in  conflict  with  one  another?  Using  the  above  information, 
the  optimal  allocation  software  was  run  and  the  results  are  printed  below. 

OPTIMAL  BUDGET  ALLOCATION  PROGRAM 


Desired  Ratio  =  7.00 


Retirement  Rate  for  Aircraft  *  0.067 
Retirement  Rate  for  Ships  =  0.033 
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40.00 

552.00 
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Number  of  Ships  vs  Time  (7  to  1  Ratio) 


It  can  be  seen  that  given  the  presented  criteria,  a  500-ship  fleet  could  not 
be  maintained. 

Now,  however,  say  it  was  decided  that  a  5  to  1  ratio  of  aircraft  to  ships 
would  be  acceptable  if  the  500-ship  fleet  could  be  maintained.  Running  the  model 
again  with  this  new  ratio  the  following  results  were  obtained.  It  should  be  noticed 
that  with  the  new  ratio  it  is  indeed  possible  to  maintain  a  500-ship  navy,  given 
the  available  budgets.  Hence,  the  model  can  be  used  to  clarify  various  inconsis¬ 
tencies  in  the  inputs  to  the  model  allowing  decision  makers  to  adjust  their 
thinking  in  such  a  way  as  to  better  achieve  their  goals . 


OPTIMAL  BUDGET  ALLOCATION  PROGRAM 


Dei  i  re<j  Ratio  =  5.00 

Retirement  Rate  for  Aircraft  =  0.067 
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#  of  Ships 


Nutober  of  Ships  vs  Time  (5  to  1  Ratio) 


SOME  CONCLUSIONS 

What  do  these  examples  say  about  the  model's  usefulness?  To  a  large  extent 
that  depends  on  the  information  the  user  has  available  to  input  into  the  model. 
If  one  has  all  the  information  given  a  few  pages  ago  as  necessary  inputs  into 
the  model  and  the  only  criteria  being  considered  is  the  force  mix  ratio,  then 
the  results  from  the  optimal  budget  allocation  program  should  give  a  very  good 
guide  as  to  budget  allocations  for  all  years  in  the  planning  period. 

But  what  if  information  about  the  future  is  unknown  (i,e.,  budgets  and  cost 
are  unknown  in  the  future)?  Need  the  methodology  presented  here  be  abandoned? 

The  answer  is  no,  if  one  is  willing  to  forecast  into  the  future  based  on  what 
has  happened  in  the  past.  This  was  the  problem  faced  when  incorporating  the 
optimal  allocation  model  into  a  larger  model  of  the  Navy.  The  methodology 
developed  to  overcome  the  problem  of  the  lack  of  future  information  is  presented 
in  the  next  section. 
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INCORPORATION  CF  THE  ALLOCATION  MODEL  IN  A  HIGHER  LEVEL  MODEL 

The  initial  purpose  of  this  project  was  to  develop  a  method  whereby  a  more 
general  model  of  naval  budgeting,  maintenance  etc.,  could  use  some  method  of 
optimal  allocation  for  its  budget  available  to  procure  aircraft  and  ships. 

Prior  to  implementation  of  this  method  the  allocation  of  the  budget  was  simply 
a  50~50  split  with  small  adjustments  for  past  orders  of  aircraft  and  ships . 

The  method  could  now  be  modified  to  entail  the  following: 

1)  Based  on  past  data  project  budget  and  cost  growth  rates 

2)  Using  this  projection  compute  estimated  budgets  and  costs  over  a 
10-year  period 

3)  Using  these  estimates  and  other  (known)  inputs,  use  the  optimal  allocation 
scheme  to  find  the  optimal  proportion  to  be  spent  on  aircraft  each  year 

in  the  planning  period. 

4)  Use  the  computed  first  year  proportion  to  allocate  the  present  year's 
budget 

5)  Repeat  the  above  process  for  each  year  in  the  model 

This  methodology  allows  a  higher  level  model  to  use  the  optimal  allocation 
scheme  without  the  need  for  future  data  (i.e.,  future  budgets  and  costs). 

Thus  a  model  which  models  only  one  year  at  a  time  with  no  lookahead  capability 
can  utilize  a  form  of  optimal  allocation. 

The  results  from  using  this  procedure  should  be  quite  good  because  previous 
years'  allocations  are  incorporated  into  the  present  year's  allocations  through 
the  backorders  of  ships  which  are  inputs  to  the  allocation  model.  Hence,  if 
many  aircraft  are  purchased  in  one  year,  this  would  be  reflected  in  the  back¬ 
orders  of  aircraft  in  the  next  year,  thus  decreasing  the  optimal  proportion  for 
this  year.  Also,  since  the  budget  and  cost  growth  rates  can  be  adjusted  each 
year,  this  method  of  allocation  should  be  able  to  provide  reasonable  results 
if  the  growth  rate  changes  are  not  too  sharp. 
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DERIVATION  OF  THE  MODEL 

At  this  point  it  would  be  prudent  to  describe  in  some  detail  the  derivation 
of  the  model  employed  to  develop  the  preceding  results.  This  section  will  first 
describe  the  derivation  of  formulas  giving  the  number  of  aircraft  and  ships  in 
any  year  given  the  proportion  of  the  budget  to  be  spent  on  aircraft  each  year 
from  the  current  year  until  the  year  in  question  and  past  data.  Next,  a  deri¬ 
vation  of  an  objective  function  for  the  model  is  presented,  along  with  its 
gradient  (which  is  used  in  the  optimization).  And  finally  a  presentation  of 
the  actual  optimization  method  is  given. 
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DERIVATION  OF  THE  NUMBER  OF  AIRCRAFT  AND  SHIPS  IN  YEAR  1 

To  begin  with  consider  the  following  definitions. 

Lets 

u  =  retirement  rate  for  aircraft  (in  decimal  form  i.e., .03  for  3%  per  yeax) 
v  =  retirement  rate  for  ships 

b^=  budget  available  for  procurement  of  aircraft  and  ships  in  year  i 
ca^=  cost  of  aircraft  in  year  i  (average) 
cs^=  cost  of  ships  in  year  i 

r^=  proportion  of  the  budget  spent  on  aircraft  in  yeax  i 
na^=  number  of  aircraft  ordered  in  year  i  (to  be  delivered  in  year  i+2) 
ns^=  number  of  ships  ordered  in  year  i  (to  be  delivered  in  year  i+5) 

A,=  number  of  aircraft  in  the  fleet  at  year  i 
S  =  number  of  ships  in  the  fleet  at  year  i 

Since  the  retirement  rate  for  aircraft  and  ships  was  assumed  to  be  constant 

and  the  leadtimes  for  aircraft  and  ships  are  2  and  5  years  respectively,  then: 

Ai  *  d-“)Ai-i  +  nai-2 

Si  *  +  nsi-5 

These  equations  say,  for  instance,  that  the  number  of  aircraft  in  year  'i  equals 

the  number  of  aircraft  on  hand  in  year  i-1  still  active  in  year  i  (lOO(l-u)^) 

plus  the  number  of  aircraft  ordered  in  year  i-2.  In  a  like  manner,  Ai  and 

Sj__1  can  be  expressed  in  terms  of  A^_^  and.  so  on  to  an<^  ^0’ 

From  this  recursive  expression  the  terms  A^  and  can  be  expressed  as 

functions  of  Aq  and  Sq  plus  the  number  of  ships  and  aircraft  ordered  prior 

to  year  i .  After  some  algebraic  manipulation  the  following  can  be  obtained: 

,  i-l 

A  =  (l-u)  A  +  (l-u)J(na  ) 

j=0 

1  i-l 

Si  =  (l-v)  SQ  +  ^  (l-v)J(ns  .  -) 

j=0  J  D 

Notice  that  the  indexes  for  na  an(i  ns-  •  c.  can  become  negative.  These 

j  ^  ~  J 

negative  indicies  represent  backorders  of  ships  and  aircraft,  ordered  before 
the  planning  period  began. 
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Now  if  these  negative  indicies  are  pulled  out  of  the  summations  the  following 
equations  are  obtained: 

Ai=(l-u)iA0+(l-u)i-1na1_1  +(l-u)i_2na0  +  ^  (l-u^na^^ 

0 

Si=(l-v)iS0+(l-v)1  1nsi_z++(l-v)i'nsi_2+(l-v)i"^nsi_2+(l-v)i_i4'nsi_1 

+  (l-v)i_5ns0  +  ^  (l-v)^ns  .  - 
j=Q  J  J 

The  terms  preceding  the  summation  in  each  of  the  above  equations  are  based 
on  actions  which  have  already  occured  and  are  in  no  way  affected  by  budget 
allocations  in  the  planning  period  (whereas  the  other  na^  and  ns^  are). 

The  terms  can  be  grouped  for  simplicity  into  functions  and  Y^  for  the 
A^  and  terms  respectively.  Hence,  with  this  simplification: 


i-3 

Ai  =  x i  +  Zi 
3=0 
i-6 

si  =  Y  f  C 
3=0 

Now,  examining  na 


nai  *  W*! 


(l-u)^na^_j_2 

(l-v)Jnsi-j-5 


and  ns^  it  can  be  seen  that: 
and  ns^  =  (l-r^b^/cs^^ 


Thus ,  the  number  of  aircraft  ordered  in  year  i  equals  the  proportion  of  the 
budget  spent  on  aircraft  (r^)  times  the  budget  (B^)  (yielding  the  amount  spent 
on  aircraft)  divided  by  the  average  unit  cost  of  aircraft  (ca^.  Likewise  for 
ships,  but  here  the  proportion  spent  on  ships  is  l-(prop.  spent  on  aircraft). 
Now,  if  these  substitutions  are  made  and  the  indexes  shifted  the  following 
are  true: 

h  ■  h  +  g  rkVcak 

Si  =  Ti  +  g  (l-v)t'k'5  Vcsk 

Hence,  if  the  b^,  ca^,  cs^#  and  Y^  are  known  the  above  equations  represent 

the  force  levels  for  aircraft  and  ships  in  terms  of  the  proportion  of  the  budget 

to  be  spent  on  aircraft  in  each  year  of  the  planning  period  (i.e.,the  r.  )  . 
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DERIVATION  OF  THE  OBJECTIVE  FUNCTION 

As  discussed  in  the  first  section  of  this  paper  the  objective  of  an 
optimal  budgeting  allocation  scheme  would  be  to  choose  proportions  of  the 
budget  to  be  spent  each  year  on  aircraft  so  as  to  keep  the  actual  ratio  of 
aircraft  to  ships  as  close  as  possible  to  some  desired  ratio.  In  other  words, 
we  want  to  find  a  set  of  proportions  (r^)  that  minimize  some  cost  function 
relating  to  the  deviation  of  the  actual  ratio  from  the  desired  ratio  in  all 
the  years  of  the  planning  period  (plus  two  years  since  the  last  year ' s  policy 
will  have  no  effect  until  then)  .  Since  deviations  above  or  below  the  desired 
ratio  can  be  considered  equally  costly  the  objective  might  be  to  minimize  the 
sum  of  the  squared  deviations  of  the  actual  ratio  from  the  desired  ratio  for 
all  years  in  the  planning  period. 

Now  Let: 

A^  =  Number  of  aircraft  in  the  fleet  at  year  i 
=  Number  of  ships  in  the  fleet  at  year  i 
k  =  Desired  ratio 

tmax  =  Number  of  years  in  the  planning  period  plus  2 

r  =  Vector  of  proportion  of  the  budget  to  be  spent  on  aircraft  (l  through 
tmax- 2) 

The  objective  then  becomes: 

tmax  A 

minimize  ^  (—5 —  -  k  ) 

i=3  bi 

Notice  that  this  summation  runs  from  i=3  to  tmax  rather  than  i=l  to  tmax.  This 
is  because  prior  to  the  third  year  of  the  planning  period  no  budget  allocation 
chosen  during  the  planning  period  will  have  any  effect  on  the  ratio  of  aircraft 
to  ships  and  hence  need  not  be  included  in  the  objective  function.  Now  let: 

A, 

^(r)  =  -  k) 

Note  also  that  the  dependency  of  f^  on  the  vector  r  is  through  the  A^  and 
terms . 
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So,  with  this  definition  the  objective  function  becomes: 

tmax  ?  tmax  _ 

minimize  (f ,(r)r  and  let  ^  (f .  (r)  r  =  D(r) 

r  i=3  1  i=3  1 

To  minimize  a  nonlinear  function  such  as  this  one  procedure  would  be  to 
choose  an  initial  starting  point  (vector  )  for  r  and  then,  since  the  opposite 
direction  of  the  gradient  of  a  function  is  the  direction  of  steepest  descent 
of  the  function,  move  along  the  opposite  direction  of  the  gradient  until  the 
desired  minimum  is  reached.  To  do  this  the  gradient  of  D(r)  with  respect  to 
the  r^  is  needed. 

DERIVATION  OF  THE  GRADIENT 

From  the  above  definition  of  D(r)  it  can  be  seen  that: 

-f§-  =  2(fi(r))4V± 

i  i=3  1  Tr" 

W 


So , 


3  i=3  Br. 

J 


If  we  represent  the  f  (r)  as  a  vector  f  and  the  \ f  1  likewise  then: 

Jr. 

=  f  •  (  ^4  . ?£tmax) 


Sr.  5r, 
J  J 


So, 


3d 

2£l  7,1 4  .  2  f tmax 

AVD  = 

**1 

3^  ar^  dr. 

•  • 

l 

•  t 

•  * 

•  • 

3D 

•  » 

. 

B  r 

_  tmax-2 

dr  3r;  ar 

tmax 


Again,  the  f's  run  from  3  to  tmax  because  no  allocation  policy  has  any  effect, 
because  of  leadtimes,  prior  to  year  3*  Also  the  r^  run  from  1  to  tmax-2  because 
proportions  cho=on  in  years  tmax-1  and  tmax  will  have  no  effect  in  the  time  frame 
minimized  over. 
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Now  the  preceding  equation  can  be  simplified  by  noticing  that  the  matrix 
is  the  transpose  of  the  Jacobian  Matrix  of  the  and  the  r±  (only  those 
applicable  to  the  problem) .  If  we  denote  this  matrix  by  J  then: 

£vd  =  jt  f 


But  what  of  the  elements  of  the  Jacobian  Matrix?  How  are  they  computed? 

DERIVATION  CF  THE  JACOBIAN  MATRIX 

A  typical  element  of  the  Jacobian  matrix  is  ^fj  . 

ari 

Now, 


111  -  ^  (  A.  .x 

Srj  S7  -k)  ■ 


(Si)' 


Further,  from  before  we  have  that: 

i-2  ... 

Ai  =  Xi  +  ^  (l-u)1_J_2b.r./ca. 

j=l  J  J  J 


1-5 


1-3-5, 


si  =  Yi  +  d  (l-v)  ■J“-?b  (l-r 
j=l  J  J 


-r  .)/cs. 


S  o  | 


^Ai  = 


for  j  > i-2 


b.(l-u)i_j~2/ca.  for  .j  <i-2 
J  J 


!fi  = 
^rj 


for  j>  i-5 


-(l-v)i"j"5b./cs.  for  j<i-5 
J  J 


combining  this  information  and  substituting  into  the  previous  equation  for  the 
partial  derivatives  of  the  ^  we  obtain  the  results  on  the  following  page. 
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ELEMENTS  OF  THE  JACOBIAN  MATRIX 


For  i=3 
Sf 


"rj 


2  = 


for  j  >  1 


(‘bi/ca1)  (l/S^)  for  j=l 


For  1=4 


bf 


4  =  < 


Sr . 
3 


\  0  for  j>2 

(l-u)4“j"V/(caj  S4)  for  j  £  2 


For  1=5 

3f, 


$rj 


0  for  j>3 

(l-u)5'j"2bj/(caj  S5)  for  j£3 


For  1  >5 
fcf 


1  = 


S1(l-u)1~j~2(b./cai)  +  Ai(l-v)1~j~5Cb./cs.) 


(si)i 


(l-u)i-j-2l,j/(c^j  Si) 


for  j  > 1-1 
for  1  £  j  £.1-5 

for  j=i-4,i~3,i-2 
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From  the  preceding  derivations  the  gradient  of  D(r)  (VD)  can  be  computed 

at  any  point.  This  computation  is  done  by  first  computing  the  Jacobian  Matrix 

as  described  then  multiplying  this  matrix  by  the  vector  of  deviations  (the  f 

T  T 

vector  described) .  This  multiplication  yields  J  f=  -§-VD  and  multiplying  J  f 
by  2  yields  VD.  From  this  point  an  optimization  algorithm  can  be  described 
to  take  an  initial  starting  point  (r  vector)  and  iteratively  improve  it  by 
adjusting  it  by  some  increment  in  the  opposite  direction  of  the  gradient  until 
the  minimum  of  D(r)  is  found, to  within  some  tolerance. 

METHOD  OF  OPTIMIZATION 


The  chosen  method  for  finding  the  minimum  of  the  D(i)  function  utilizes 
a  very  simplified  gradient  seach  method.  The  algorithm  used  can  be  described 
as  follows. 

1)  Set  an  initial  proportion  vector  r  (Chosen  to  be  all  .5) 

2)  Set  an  initial  increment  to  adjust  the  r^  elements  (Chosen  to  be  .4) 

3)  Compute  the  "Jacobian  Matrix  sis  described  in  the  previous  section  (j) 

4)  Compute  the  vector  of  deviations  (f)  given  the  current  r  vector 

5)  Compute  J  f  (this  is  j  the  gradient  of  the  objective  function) 

6)  Since  the  objective  is  to  minimize  the  objective  function  we  would 
like  to  move  in  the  opposite  direction  of  the  gradient.  So  the 
algorithm  implemented  does  the  following. 

If  m  , 

(j  f )  <  0  then  let  r^  =  r^  +  increment  (if>l  r^=l) 

(j^f)  >  0  then  let  r^  =  r^  -  increment  (if<0  ri=0) 

(j^f^  =  0  then  let  r^  =  r^ 

7)  Now,  check  if  all  the  r,  are  equal  to  either  the  preceding  r^  or  the 
r,  preceding  that.  If  they  all  are,  then  the  solution  is  oscillating 
between  two  points  and  so  the  optimum  should  be  somewhere  in  between 

the  two  vectors  that  are  being  oscillated  between.  So,  if  this  condition 
is  true  the  increment  is  cut  in  half  and  the  algorithm  proceeds. 

8)  The  process  described  above  continues  until  the  increment  is  smaller  than 
some  set  tolerance  level.  (Chosen  to  be  .001.)* 


* 

The  .001  tolerance  level  was  chosen  because  it  was  felt  that  changes  to  the 
r  vector  smaller  than  that  were  not  really  significant  enough  to  be  taken  into 
account.  This  increment  can  be  easily  changed  if  it  were  desireable  to  do  so. 
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SOFTWARE  FOR  OPTIMAL  ALLOCATION 

This  section  presents  the  software  developed  for  the  optimal  budget 
allocation  problem  discussed  in  this  paper.  The  software  was  written  in 
both  standard  FORTRAN  and  FORTRAN-77  for  implementation  on  an  IBM  and  a 
VAX  II/78O  system  respectively.  The  software  for  both  systems  has  been 
validated  but  only  the  software  for  the  IBM  implementation  is  presented 
here . 
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THE  DEVELOPED  SOFTWARE 

The  software  developed  to  do  the  optimal  budget  allocation  is  written  in 

standard  FORTRAN  and  consists  of  4  primary  routines  and  a  routine  written  as 

a  linkage  between  the  optimal  allocation  routines  and  another  model.  A 

description  of  the  routines  follows. 

ALLOCA  -  This  routine  performs  the  optimization  described  in  the  previous  section 
of  this  paper.  It  should  be  noted  that  this  routine  may  be  replaced  by 
other  (possibly  more  sophisticated)  nonlinear  optimization  routines  if 
the  convergence  of  the  outlined  methodology  is  considered  too  slow.  I 
don't  consider  its  convergence  is  bad,  but  bad  is  a  somewhat  relative  term. 

JACOBI  -  This  routine  computes  the  Jacobian  Matrix  (J)  of  the  problem  using 
the  equations  developed  in  the  preceding  section  of  this  paper. 

JTRANF  -  This  routine  multiplies  the  vector  of  deviations  of  the  actual  ratio 
from  the  desired  ratio  by  the  transpose  of  the  Jacobian  Matrix  computed  by 
JACOBI.  This  yields  the  vector  JTf,  \  the  gradient  of  the  objective  function. 

C0M5TR  -  (stands  for  compute  strength)  This  routine  computes  the  number  of  air¬ 
craft  and  ships  in  any  year  of  the  planning  period  plus  two  years ,  given 
past  data  and  a  vector  of  proportions  of  the  budget  to  be  spent  on  aircraft 
in  each  year  of  the  planning  period. 

BUDALL  -  This  is  a  linkage  routine  between  a  higher  level  model  and  the  above 
optimal  allocation  routines .  Its  primary  function  is  to  simplify  the 
inputs  to  the  model — requiring  the  calling  program  to  specify  only  budget 
and  cost  growth  rates,  and  to  place  the  calling  programs  inputs  into  a  block 
of  common  storage  accessed  by  all  the  optimization  routines.  This  routine 
should  be  easily  modified  to  suit  the  users  precise  specifications.  For 
instance,  if  growth  rates  are  not  constant,  the  calling  program  could  pass 
vectors  containing  the  budgets  and  costs  for  all  years  and  then  have  BUDALL 
put  these  inputs  into  the  appropriate  location  in  COMMON  storage 
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CALLING  INSTRUCTIONS 

As  presently  formulated  the  calling  sequence  for  the  optimization  routines 
is  via  the  routine  BUDALL.  To  access  the  routines  the  following  FORTRAN  state¬ 
ment  is  needed: 


CALL  BUDALL  (CURAIR .CURSHI, BUDNOW.CANOW.CSNOW, 

1  BUDGl.CAGl ,CSG1  ,BACKA, BACKS ,TIME,X,Y, 

1  RATCON.R) 

Where  the  l’s  represent  FCRTRAN  continuation  cards  and  the  parameters  are 
defined  as  follows  * 

CURAIR  -  The  number  of  aircraft  currently  in  the  fleet  (in  year  zero) , INTEGER 

CURSHI  -  The  number  of  ships  currently  in  the  fleet  (in  year  zero) , INTEGER 

BUDNOW  -  The  budget  available  for  procuring  aircraft  and  ships  in  year  1 

CANOW  -  The  average  unit  cost  of  aircraft  in  year  1  of  the  planning  period 

CSNOW  -  The  average  unit  cost  of  ships  in  year  1  of  the  planning  period 

BUDG1  -  The  growth  of  the  available  budget  each  year(i.e.,3^  growth,budgl=,03) 

CAG1  -  The  growth  rate  for  aircraft  cost  (similar  formulation  as  for  budgl) 

CSG1  -  The  growth  rate  for  ship  cost  (same  as  for  budgl) 

BACKA  -  Vector  of  length  2  giving  the  backorders  of  aircraft .INTEGER 

with  those  ordered  1  year  ago  in  BACKA(l) 
with  those  ordered  2  years  ago  in  BACKA (2) 

BACKS  -  Vector  of  length  5  giving  the  backorders  of  ships , INTEGER 
with  those  ordered  1  year  ago  in  BACKS (l) 

with  those  ordered  5  years  ago  in  BACKS (3) 

The  number  of  years  in  the  planning  period, INTEGER 
The  retirement  rate  for  aircraft  (3^  retirement  rate,  x=.03) 

The  retirement  rate  for  ships  (same  as  above) 

The  desired  ratio  of  aircraft  to  ships 

A  vector  of  length  30  in  which  the  vector  of  optimal  proportions 
will  be  returned,  where  R(i)=  the  proportion  of  the  budget  to  be 
spent  on  aircraft  in  year  i 

If  it  is  desired  to  modify  BUDALL  the  necessary  modification  to  the  above  call 

statement  will  need  to  be  made . 


TIME  - 

X 

Y 

RATCON  - 

R  - 
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COMMON  BLOCK 

Regardless  of  the  modification  made  to  the  routine  BUDALL  to  suit  the 
users  specifications ,  BUDALL  must  place  values  of  the  input  parameters  into 
the  COMMON  block  ALLOC  within  that  subroutine  (this  is  a  labelled  C OMMON 
block  in  FORTRAN).  The  following  parameters  must  be  put  in  common.* 


ITIME  -  The  number  of  years  in  the  planning  period,  INTEGER 

B  -  A  vector  of  length  30  containing  the  available  budget  for  each 

year  in  the  planning  period 

CA  -  A  vector  of  length  30  containing  the  unit  aircraft  cost  for  each 
year  in  the  planning  period 

CS  -  A  vector  of  length  30  containing  the  unit  ship  cost  for  each 
year  in  the  planning  period 

P  -  ( 1-Retirement  rate  for  aircraft) 

Q  -  (l-Retirement  rate  for  ships) 


BACKA  -  Same  definition  as  previously  given, INTEGER 
BACKS  -  "  "  "  INTEGER 
CURAIR-  "  "  "  INTEGER 
CURSHI-  "  ’’  "  INTEGER 
RATCON-  "  " 

TOL  -  The  tolerance  acceptable  for  minimization  (either  .001  or  .002 
should  be  acceptable) 


The  common  statement  should  take  the  following  form. 


COMMON/ALL0C/lTIME,B(30),CA(30),CS(30) ,P,Q,BACKA(2),BACKS(5) 
1  .CURAIR, CURSHI ,RATC0N, TOL 


The  above  variables  also  represent  most  of  the  important  variables  for  each  of 
the  routines.  For  further  information  on  the  software,  one  should  either  examine 
the  FORTRAN  code  and  the  comments  contained  therein,  or  refer  to  the  flow  diagrams 
on  the  following  pages. 


It  should  be  noted  that  in  the  code  for  BUDALL  there  is  a  slight  modification 
to  some  of  the  variable  names  in  the  common  block.  This  is  because  it  was 
desired  to  use  the  same  variable  names  in  the  parameter  list  passed  as  were 
used  throughout  the  program.  FORTRAN  doesn't  allow  a  variable  to  be  in  both 
a  parameter  list  and  a  common  statement,  so  to  circumvent  this  problem  the 
variable  in  this  routine  alone  had  their  names  changed  from  what  is  above. 
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HIERARCHY  DIAGRAM  FOR  OPTIMAL  BUDGET  ALLOCATION 

PROGRAM 
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Significant  Variables:  (Not  Explained  Elsewhere) 

TMAX  -  #  of  Years  in  the  Planning  Period  plus  2  Years,  INTEGER 
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Significant  Variables:  (Not  Explained  Elsewhere) 

r(i)  -  Vector  of  length  JO  to  contain  the  proportion  of  the  budget 
to  be  spent  each  year  on  aircraft 
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Significant  Variables;  (Not  Explained  Elsewhere) 

air(i)  -  vector  of  length  30  to  contain  the  #  of  aircraft  in  the  fleet  at 
year  i  in  air(i) 

ships (i)  -  vector  of  length  30  to  contain  the  #  of  ships  in  the  fleet  at 
year  i  in  ships (i) 
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Significant  Variables:  (Not  Explained  Elsewhere) 

a(i)  ~  a  vector  of  length  30  to  contain  the  #  of  aircraft  in  each  year 

s(i)  -  "  ”  "  ships 

JTF(i)  -  vector  of  length  30  to  contain  the  elements  of  the  vector  J^F 
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Significant  Variables:  (Not  Explained  Elsewhere) 


a(i) 

s(i) 


-  vector  of  length  30  to  contain  the  #  of  aircraft  in  the  fleet  at 
year  i  in  a(i) 

-  vector  of  length  30  to  contain  the  #  of  ships  in  the  fleet  at 
year  i  in  s(i) 
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CONCLUSIONS 

Given  the  preceding  results,  what  conclusions  can  be  drawn?  The  model 
of  the  budget  allocation  process  as  developed  and  solved  allows  for  the  op¬ 
timal  allocation  of  a  budget,  over  some  time  span,  toward  procurement  of  air¬ 
craft  and  ships  so  as  to  keep  the  ratio  of  aircraft  to  ships  (force  mix)  as 
close  as  possible  to  some  desired  level,  regardless  of  what  that  level  may  be. 
The  developed  software  does  this  by  minimizing  the  sum  of  the  squared  de¬ 
viation  of  the  actual  ratio  from  the  desired  ratio  for  all  years  in  the  planning 
period . 

The  results  presented  in  this  study  are  generally  quite  problem  specific, 
dealing  with  the  particular  problem  of  allocating  a  budget  between  purchases 
of  aircraft  and  purchases  of  ships  so  as  to  hold  the  force  mix  constant.  If 
one  is  seeking  answers  to  such  questions  or  if  one  is  interested  in  studying 
the  interactions  of  the  model’s  inputs  or  their  effect  on  force  levels,  the 
model  and  solution  method  presented  in  this  paper  should  serve  as  a  useful  tool. 
If,  on  the  other  hand,  one  is  looking  for  a  methodology  for  solving  other,  sim¬ 
ilar,  allocation  problems,  the  approach  presented  here  may  be  useful  as  a  guide, 
because  of  its  simplicity. 

Throughout  this  study,  it  was  felt  that  one  prime  objective  in  the  develop¬ 
ment  of  the  foregoing  model  was  that  the  model  should  remain  as  uncomplicated 
as  possible  while  still  yielding  a  correct  solution.  This,  to  a  large  extent, 
was  accomplished. 

Likewise,  the  prime  objective  of  the  study:  The  optimal  allocation  of  a 
budget  toward  procurement  of  aircraft  and  ships  in  such  a  way  as  to  keep  the 
force  mix  constant,  was  achieved. 
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Appendix  I  -  Some  further  results 

On  the  following  two  pages  are  the  results  of  two  runs  of  the  optimal  budget 
allocation  software.  The  inputs  to  these  runs  are  the  same  as  for  the  first 
example  in  section  one  except  that  the  desired  ratio  for  these  runs  is  8  for 
the  first  and  5  for  "the  second.  These  results  are  presented  to  further  dem¬ 
onstration  of  the  ability  of  the  optimal  allocation  software  to  achieve  its 
objective  of  keeping  the  force  mix  between  aircraft  and  ships  at  the  desired 
level. 
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Appendix  II  -  FORTRAN  code  for  optimal  allocation 

Routines  listed: 

BUDALL  -  Linkage  Routine 

ALLOCA  -  Optimal  Allocation  Optimization  Routine 

JACOBI  -  Computes  the  Jacobian  Matrix  of  the  Problem 

T 

JTRANF  -  Computes  J  f  (the  gradient  of  the  objective  function) 

C0M5TR  -  Computes  the  number  of  aircraft  and  ships  in  the  fleet  at  any 
year  given  a  particular  allocation  policy  over  some  planning 
period 
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C  +  *'*  +  *’*+***'*+  ’*  +  **'*  +  T  **+'*  +  +  +**+**'*'**+'***  m  n*  +*+*++*+  +  +***+  *  + 

c  * 

C  *  ROUTINE  3UDALL  -  THIS  IS  A  LINKAGE  ROUTIN'  E  BETWEEN  * 

-  *  liiE  OPTIMAL  BUDGET  ALLOCATION  ROUTINES  AND  THE  * 

0  *  USER'S  SOFTS*  ARE .  * 


C 

r 


\  ^***^**W  *+★***  +  It*  ~  ar***-*-**-****,,  ♦  »★*-«■*»**-»-*'*•+*****-**  *-**■#* 

SUBROUTINE  BUDALL  (  C  UP  A  I  R  ,  C  UR  S  H 1 , 3  U  D  NO  W  ,  C  A  li  0  W ,  C  S  M  C  W , 

1  EUDG 1 ,  C  AG  1 ,0301 , BACK A, BACKS , T IKE , X , Y , 

1  EATCCM,R) 

INTEGER  BACK  A ( 2 ) , 8 AC  NS ( 5 ) , C  UR A  I K,CORSH I , TM  A X ,TI ME 
INTEGER  CUBS  5I ,  CURP.SB,  BA,  BS 
REAL  ROC) 

CO  MUON/ ALLOC/ITIM  E  ,  B  (3 C  )  , C A  (  3 0  )  ,CS  (  3 0  )  ,  P  ,0 , 3  A (  2  )  ,  BS  (  5  )  , 

1  CUES  AI , CU  3  3S  H , R A? , TOL 

?  AND  Q  EQUAL  THE  PROPORTION  OF  THE  FLEET  ACTIVE  IN  YEAS  I 
STILL  ACTIVE  IN  YEAR  1+1  ( 1 -? ET IE E M £ NT  RATE) 

P=1-'( 

q=i-  y 


THIS  IS  THE  TOLERANCE  LEVEL  FOR  THE  OPTIMIZATION  ALGORITHM 


TO L=. 0020 
3UDGBW=1 +3UDG 1 
CAGB W= 1 +CAG1 
C3GR  W=1+CSG1 


C  INFORMATION  FROM  CALLING  PROGRAM  PUT  INTO  COMMON  STORAGE 
r 

CU3R  AI=CUR AI c 
CUSRSH=CUHSHI 
RAT= AATCON 
3C  1  )  =  B  U  0  N  0  W 
C  A ( 1 ) =C A NO Vf 
C3( 1 ) =CSNOW 
r 

c  set  tmax  -  a  of  years  over  which  the  objective  is  minimized 
tiak=tim::-+2 


C  FILL  IN  BUDGET  AND  COST  VECTORS  FOR  ALL  YEARS  IN  THE  PLANNING  PERIOD 

DO  100  1  =  2,  TMAX 

3 ( i )=rudgrw+ s( r-i ) 

C  A ( I ) =C AGF V’C\(I-1 ) 

C3(I)=CSGRW'*CS(I-1) 

131  GONTIVUE 

r 

C  1  OR  £  PARAMETERS  PUT  INTO  COMMON 

r* 

3A( 1  )  =  3  A  C  X  A  ( 1  ) 

3  A  (  2  )  =  3  A  C  X  A  (  2) 

DO  200  1=1,5 

3  3  ( I )  =  3  ACK.F  ( I ) 

230  CONTINUE 

I T I M  F = T 1 1 E 
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ALL  OPTIMIZATION  ROUTINE  ALLOCA  TO  COMPUTE  OPTIMAL  PROPORTION  OF 
HE  OJDCET  TO  3E  SPENT  OH  AIRCRAFT  FOP  EACH  TEAR  IN  THE  PLANNING 
C  PF3IQ0.  THE  RESULT  IS  RETURNED  IN  THE  VECTOR  R  WHERE  THE  OPTIMAL 
C  PROPORTION  FOR  YEAR  I  IS  RETURNED  IN  *(I). 

5*  CALL  ALLOCA  ( R , TMA  X ) 

59  RETURN 

50  END 
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*-<•-*•****-**  *  *  if  w  ft  +  +  ’ 


C 

n 

r 


ROUTINE  ALLOC 
ATI 0."  OF 
AIRCRAFT 
BEHIND  T 

THE  ARRAYS  J( 
THE  JACOBIAN 


*-*★*  ***•★•*  +  ***  a****^******-**** 

* 

A  -  THIS  ROUTINE  CARRIES  OUT  THE  QPTIMIZ-  * 
THE  PROPORTION  0?  THE  BUDGET  SPENT  ON  * 

EACH  YEAR.  FOR  DETAILS  ON  THE  THEORY  * 

HE  A.  LG  OR  IT  H  M  SEE  THE  ACCOMPANYING  PAPES 

* 

30,3 0)  AND  JTF ( 30 )  ARE  USED  TO  CONTAIN  THE  * 
TAT SIX  AND  GRADIENT  (J’F)  RESPECTIVELY  * 


IT  SHOULD  BE 
OTHER  NON- LINT 
IZING  THE  ROD- 
GRADIENT  OF  T! 


NOTED  THAT  THIS  ROUTINE  MAY  BE  REPLACED  3Y  * 
EAR  OPTIMIZATION  ROUTINES  WHILE  STILL  UTIL-* 

TINES  JACOBI  AND  JT3ANF  TO  COMPUTE  THE  * 

NON-LINEAR  OBJECTIVE  FUNCTION.  * 

*  •* 

**-»**4**»  »*  +  -*■»•»***•***»*■*•**  »*■».**-..**  *■»  ***■»*•-**•-*■*■  -* 


6  1 
62 
6  3 
64 

66 
6  6 


SUBROUTINE  ALLCCA  ( R  ,TMAX ) 

REAL  J(30,30),JTF(30),P,Q,IHC,a(30),R1(30),B2(30) 

IN  IE  GSR  E  ACX  A, BACKS, TM  AX ,C USAI S ,CURS  HI 

CO  MMON/ALLOC/ITIME,  3  (30  )  ,CA(3C)  ,CS  (  30  ,  P  ,Q  , B  ACK  A  {  2 )  ,  3  AC KS(  5  ) 
1  ,CURATR,C;J53HI,RATC0N,T0L 

LOGICAL  FLAG 
INTEGER  COUNT 


INITIALIZE  PROPORTION  TC 


67 

DO  100  1  = 

6  8 

5 ( I )=L 

6  Q 

3 1  (  I )  = 

70 

R  2  ( I )  = 

71 

100  CONTINUE 

72 

INC=.4 

73 

FLAG-. TRU 

7a 
75 
75 
77 
7  c 


75 

80 

3  1 


C  LOOP  UNTIL  THE  INCREMENT  THAT  THE  R(I)  ARE  ADJUSTED  BY  IS  LESS 
C  THAN  I  HE  TOLERANCE  (THEN  FLAG  SET  FALSE) 

20 0  IF  (.NOT. FLAG)  GO  TC  1000 
DO  25C  N« 1 , ITT  ME 
R2(N)=R 1(N) 

R 1  (  N  )  =  F  (  N  ) 

250  CONTINUE 

Q 

C  CALL  ROUTINES  JACCfil  AND  JTRANF  TO  COMPUTE  THE  GRADIENT 

r* 

CALL  JACOBI  (J, P  ,  T  M  A  X  ) 

CALL  JTRANF  (  J  ,  JTF  ,  R  ,  Tff  A  X  ) 

COU  M?=0 


•  92 
8? 
H  a 
.  8  5 
36 


C  COMPUTE  A  NEW  ?  VECTOR  AND  TEST  IF  SOLUTION  OSCILLATING.  THIS  IS 
C  DONE  BY  TESTING  WHETHER  THE  ELEMENTS  OF  THE  CURRENT  SOLUTION  EQUAL 
C  EITHER  THE  PREVIOUS  SOLUTION  OR  THE  ONE  PRECEDING  THAT. 

C 

DO  300  M= 1 , ITI ME 

IF  (JTF(N).LT.O.O)  R< N ) = R ( N) + I NC 
IF  ( JTF(N).GT.r.O)  R(N)=R(N)-INC 
IF  ( P(X) .LT.O.G)  R ( N ) =0 . 0 
IF  (R(N).GT. 1.0)  PCS )=1.0 
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87 

DI  F 1=ABS  (  8  (  \'  )  -P.1  (  S  )  ) 

8  8 

DIF2=ARS(P(N)-B2( N) ) 

89 

IF  ( (DIF1.LT. 0.000001 ). CP. 

1  COUNT=CO(JKT  +  1 

90 

30  0 

CONTI  NOE 

9  1 

IF  (COUNT. E.Q.ITI2E)  TNC=INC/2 

92 

IF  (INC.LT.TOL)  FLAG=. FALSE. 

9  3 

CO  TO  200 

94 

1000 

CONTI NUE 

95 

EETUS'I 

96 

END 

.0  00001  )  ) 
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C  ****★•»»»-»■»«■+'»***■»»*«(•»***♦-»*★•****■****-»«»****'»*****■»**■**,-».+  *** 

C  *  * 

C  ‘  ROUTINE  JACOBI  -  THIS  ROUTINE  COMPUTES  THE  JACOBIAN  * 
C  *  MATRIX  OF  THE  PEOPLE!?  AS  DESCRIBED  I N  THE  AC  COMP-  * 

C  *  AN’YIMG  PAPER.  * 

C  *  + 

C  *  THE  ROUTINE  IS  PASSED  A  VECTOR  8  CONTAINING  THE  CURRENT* 
C  *  SOLUTION  (VECTOR  OF  OPTIMAL  PROPORTIONS  * 

C  *  * 

C  ***■********'*■•»***•'*****'**»■♦*■*'»«*  +  ***-*»*■*******»**■»******•**-»■* 


91  SUBROUTINE  JACOB I ( J , 3 ,TH A X ) 

98  INTEGER  3 AC KA , BACK S , C UR AIE, CU8SH I 

QQ  INTEGER  A IR ( 3C ) ,SHIPS(30) ,  TM  AX 

10p  REAL  J(  30,30,8(30  ) 

1C  1  COMMON/ ALLOC/ITIME, 3 ( 30 ) ,C A ( 30 ) , CS( 3 0 ) , P , 0 , BACK A ( 2 )  , BACKS (5) 

1  ,  CUR  A  IE,  CUR  SHI,  RATC  C.M  ,  TO  L 

C 

C  CALL  ROUTINE  COMSTE  TO  DETERMINE  THE  NUMBER  OF  AIRCRAFT  AND  SHIPS 
C  IN  THE  FLEET  AT  ANY  YEAR  IN  THE  PLANNING  PERIOD  GIVEN  THE  CURRENT 
C  SOLUTION 
r 

102  CALL  COSSTB  ( AIR , SHIPS , 8 , TSA X) 


ZERO  OUT  THE  J  MATRIX 


10  3 

r 

DO  20  1=1,30 

,104 

DO  10  N =1,30, 5 

10  5 

J( I,N )=C.O 

10* 

J(I,N+1)=0.0 

.107 

J(I , N+2 )=0 . 0 

103 

J(I,N+3)=0.0 

109 

J(I  ,N  +  4)=0.0 

11C 

10 

CONTINUE 

1  1  1 

20 

CONTINUE 

c 

C  COMPUTATION  OF  THE  JACOBIAN  MATRIX  ELEMENTS 


112 
113 
1  1  4 
11* 
116 

117 

118 

119 

120 
121 
122 

123 

124 

125 

126 
127 
12° 
125 

130 

131 

132 


J(3,1 )=3( 1)/(CA(1)*SHIPS(3>) 

J  (4 ,1  )  =  (?'3(  1  )  )/(  CA  (  1  )*3TJI?S  (4  )  ) 

J(5,1)=(?‘P*B(1))/(CA(1)‘SHIPS(5)) 

J(3,2)=0.0 

J(4,2 )  =  8<2)/(CA(2) ‘SHIPS (4  )  ) 
J(5,2)=(P*3(2) ) / ( C A ( 2 )‘ SHIPS ( 5 ) ) 
J(3,3)=0.0 
J(4,3 )=0. 0 

J(5,3 )  =  8( 3)/(CA( 3) ‘SHIPS ( 5) ) 

DO  800  1=4  ,  ITI  xME 
J( 3,1) =0.0 
J( 4, I)=0.0 
J( 5, I) =0.0 
800  CONTINUE 

DO  1000  1=6 , TM AX 
NUM=I-5 

DO  900  X  =  1  , N U M 

A=FLOAT( A I 8 ( I ) ) 
S=FLOAT(SHIPS(I.)) 

X  =  P“(I-X-2) 

Y  =  Q“  ( I-K-5) 


41  - 


APPENDIX  II  -  Continued  T-475 


133 

J(I,K)=((  (S*X*3(K)/Cr.  (TO  )  +  (A*T*B('0/CS(X>  )  )/  (S*S)  ) 

134 

900 

CONTINUE 

135 

NUN 1=T- 4 

-136 

NUN  2  =  1- 2 

137 

DO  95u  K=NU* 1 , NUH2 

138 

J(I,X)  =  ((  P**(I-K-2)  )  *BOO  )/(CA(X)*SHIPS(I  )  ) 

.  139 

950 

CONTINUE 

140 

1 0  0  0 
c 

C  HER 

CONTINUE 

S  THE  MATRIX  IS  SHIFTED  SO  THAT  THE  PARTIAL  DERIVATIVE  OF 

C  THE  3  YEAR'S  DEVIATION  WITH  RESPECT  TO  THE  FIRST  TEAR'S  PRC  E- 
C  PORTION  IS  FOUND  IN  THIS  AIDS  IN  THE  CONFUTATION  Of 

C  JTF(I)  THE  GRADIENT  OF  THE  OBJECTIVE  FUNCTION 


r 


141 

DO  1500  1=1,28 

142 

DO  1250  K=1 ,30 

143 

J(I,X)=J( I+2,K ) 

144 

1250 

CONTINUE 

145 

150  0 

CONTINUE 

146 

2000 

CONTINUE 

147 

RETURN 

14  8 

END 
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t"  *  .* 


c 

r* 

C 

c 

c 


*  ROUTINE  JTRANF  -  THIS  ROUTINE  TAKES  THE  TRANSPOSE  CE  * 

*  THE  J  MATRIX  COMPUTED  BY  THE  ROUTINE  JACOBI  AND 

’  POST- MULTIPLIES  IT  BY  THE  VECTOR  OF  DEVIATION  OF  * 

*  THE  ACTUAL  RATIO  FROM  THE  DESIRED  RATIO  GIVEN  THE  * 

*  CURRENT  SOLUTION.  THIS  YIELDS  THE  GRADIENT  OF  THE  * 
OF  THE  OBJECTIVE  FUNCTION. 

*  * 
•**-**■-**■*-**  +  *-*-*++■**•*  +  *■*■-*★  -w-t r  +  +*  +  +  +-*+■  *ie-k'irle'w*-k  +  wir*+-k-Kir  +  *■■*+**  -k  *  * 


149 

SUBROUTINE  JTRANF  ( J  ,  JT F  ,  R , TB AX  ) 

ISO 

INTEGER  BACKA,B ACXS, C UR A I E, CURS  HI, T MAX 

IS  1 

INTEGER  A(3C),S(3G) 

152 

REAL  J(3C,30),JTF(30),K(30) 

1c-  3 

CO  Mil  ON  /  ALLOC  /  IT  IN  E  ,  B  (  30  )  ,CA  (  30  )  ,CS  (30),  P  ,0 ,8  ACKA(  2  )  ,  BACK  S(  5  ) 

1 

,CUSAI E,CUPSHI ,PATCON ,TOL 

154 

CALL  COM ST S  ( A ,S , P ,T R A X ) 

155 

DO  5  00  1= 1 , IT I M E 

156 

JTF (I  )  =  0.C 

157 

DO  250  K=  1  ,  IT 1 2 E 

158 

A 1  =  FLO  AT(A(K  +  2)> 

159 

S1=FL0AT(S( K+2 ) ) 

160 

JTF( I)=JTF( T ) + (J ( K ,1 ) '( ( A1/S 1 ) -RAT CON ) ) 

161 

25  0 

CONTINUE 

162 

5  3  0 

CONTINUE 

16  3 

10  3  0 

CONTINUE 

154 

RETURN 

16  5 

EN  0 
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16  6 

167 

168 
169 


r  ■* 

C  * 

C  x 
C  * 
C  * 

c  * 
c  * 
c  * 
c  * 
c  * 
c  ■* 
c  * 
c  * 
c  * 
c  * 
c  *  - 

r 


BOUTINS  COHSTR  -  THIS  ROUTINE  COMPUTES  THE  NUMBER  OF 
AIRCRAFT  AND  SHIPS  IN  THE  FLEET  FOR  ANY  YEAR  IN 
THE  PLANNING  PERIOD  PLUS  TWO  YEARS  GIVEN  THE  FOL¬ 
LOWING  INFORMATION: 

CURRENT  NUMBER  OF  AIRCRAFT  AND  SHIPS 
BACKORDERS  OF  AIRCRAFT  1  AND  2  YEARS  AGO 
BACKORDERS  OF  SHIPS  1,2, 3, 4  AND  5  YEARS  AGO 
BUDGETS  FOR  PROCUREMENT  FOR  ALL  YEARS  IN  THE 
PLANNING  PERIOD 

AVERAGE  UNIT  COSTS  OF  AIRCRAFT  AND  SHIPS  FOR 
ALL  YEARS  IN  THE  PLANNING  PERIOD 
\  VECTOR  R  CONTAIN  THE  PROPORTION  OF  THE  BUD¬ 
GET  TO  BE  SPENT  ON  AIRCRAFT  IN  EACH  TEAR 


it  *  tr 


*  *  *  *  VT  +  *  *  ir  *  *  •*  ir 


SUBROUTINE  COKSTS  (A,S,R,T3AX) 

INTEGER  3ACKA,BAC;<S,CURAIR,CUSSHI,A  (  30 ) ,S( 30 ) ,TM A  X 
REAL  E ( 3 0  ) 

COMB  ON / A  LLOC/ITIM  E ,3(30),CA(30) , CS ( 3 0 ) ,P ,Q , 3 ACX A ( 2 ) , 0  AC K S ( 5 ) 

1  ,CURAIR,CURSHI,RATCCN,TCL 

FOR  FURTHER  INFORMATION  ON  THESE  COMPUTATIONS  SEE  ACCOMPANYING  PAPER, 


17C 

A ( 1 )  a I  FIX  (  ( P*  CUB  AI R  )  +3  ACSCA  (  2)  ) 

171 

S( 1)=IFIX( (Q*CURSHI)+BACK3(5) ) 

1^2 

A ( 2  )  =1 F I X  ( ( P*  A  ( 1 ) )+BACKA( 1  )  ) 

173 

S(2)-IFIX( ( Q  *  S  (  1  ) ) + BACKS ( 4  )  ) 

174 

S(  3)  =1  FIX (  (Q'SUJ  )  +•  B  A.C  K  S  (  3 )  ) 

175 

S(  4)=IFIX(  (Q'-SC  3)  )+3ACKS(2)  ) 

176 

3(5) =IFIX( (Q'S( U) )+BACKS( 1 )  ) 

177 

DO  500  1=  3 ,T  v  ?  X 

178 

A( I)=IFIX{ (P*A(I-1 ) ) +(3(1-2 ) 

17a 

50  0 

CONTINUE 

180 

DO  800  I=6,T«AX 

181 

S(I)=IEIX((Q*S(I-1) )+((1-R(I 

1«2 

90  0 

CONTINUE 

1  S  3 

1 0  0  0 

CONTINUE 

184 

RETURN 

18  5 

END 
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